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We derive the fuU counting statistics of charge transfer through a voltage biased superconducting 
junction. We find that for measurement times much longer than the inverse Josephson frequency, 
the counting statistics describes a correlated transfer of quanta of multiple electron charges, each 
quantum associated with the transfer of a single quasiparticle. An expression for the the counting 
statistics in terms of the quasiparticle scattering amplitudes is derived. 
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Due to the discrete nature of electric charge, the cur- 
rent in mesoscopic conductors generally fluctuates. Over 
the last decade, there has been an increasing interest, 
theoretical as well as experimental fj, in the physics of 
current fluctuations. Most studies have been focused on 
noise, the second moment of the fluctuations, but re- 
cently a considerable interest has been shown for the 
full distribution of charge fluctuations, the full counting 
statistics (FCS) % A variety of theoretical approaches 
to the FCS, ranging from quantum mechanical [^Q, via 
quasiclassical U to classical 0, have been developed. 
The third moment of current fluctuations was very re- 
cently measured 0, opening up the road, as well, to ex- 
perimental investigation of the higher moments of the 
fluctuations. 

Noninteracting electrons in purely normal conduc- 
tors are transferred one by one In normal- 
superconducting junctions, the charge transfer mecha- 
nism across the normal-superconducting interface, at en- 
ergies below the superconducting gap, is Andreev reflec- 
tion. As a consequence, the FCS include terms describ- 
ing correlated transfer of pairs of electrons Q . Recently, 
Belzig and Nazarov studied the FCS in superconduct- 
ing junctions with a flxed phase difference between the 
superconducting electrodes. They found that the classi- 
cal interpretation of the FCS, the probability to transfer 
a given number of electrons across the junction during the 
measurement, could imply negative probabilities. Cou- 
pling the junction to a detector 0, they showed that 
this resulted from an attempt to interpret the phenom- 
ena of supercurrent with classical means. 

In voltage biased superconducting junctions, the phys- 
ical situation is quite different. Due to the applied volt- 
age bias, the superconducting phase difference oscillates 
with the Josephson frequency 2eV/h, giving rise to both 
dc and ac-components of the current. For measurement 
times much longer than the inverse Josephson frequency. 



only the dc-current, which is dissipative, contributes to 
the net charge-transport. Microscopically, the charge is 
transported between the two superconductors via coher- 
ent multiple Andreev reflections (MAR) 11]. The cur- 
rent has been studied in various junctions, both theo- 
retically ^3 and experimentally ,13J. Recently, also the 
noise was studied [14, 1^ J£J • 

In this paper we present the FCS of charge transfer 
through a voltage biased junction, in terms of the am- 
plitudes for quasiparticle scattering. Each quasiparticle 
scattering process results in an integer number of electron 
charges being transferred across the junction. As a con- 
sequence the FCS can be interpreted in classical terms. 
At temperatures much lower than the superconducting 
energy gap A, many-quasiparticle scattering processes 
are exponentially suppressed, resulting in a simple prob- 
ability distribution, containing only the probabilities for 
single quasiparticle scattering. This distribution repro- 
duces known results for dc-current and zero-frequency 
noise. We discuss in detail the third cumulant for single 
channel junctions and diffusive junctions shorter than the 
superconducting coherence length. 

We consider a superconducting junction consisting of 
two superconducting reservoirs connected via a normal, 
mesoscopic conductor (see Fig. For simplicity of no- 
tation, we consider a junction with a single transport 
mode, the multi-mode generalization is discussed below. 
A voltage V is applied between the two reservoirs. 

The single-particle wavefunctions in the junction, solu- 
tions to the time dependent Bogoliubov-de Gennes equa- 
tion, are scattering states labelled by the incoming quasi- 
particle type, injection energy and reservoir. The scatter- 
ing states are superpositions of amplitudes for quasiparti- 
cles at energies ±neV from the injection energy, counted 
from the local chemical potential in each contact. The 
amplitude for an incoming quasiparticle of type a at en- 
ergy Em to exit the junction as a quasiparticle of type /3 
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FIG. 1: Left: A schematic picture of the junction. Right: 
Multiple Andreev reflection processes in the voltage biased 
superconducting junction. An electron-like quasiparticle in- 
jected at an energy Eo can undergo multiple Andreev reflec- 
tions before being emitted as an a-type (e or h) quasiparticle 
at an energy En — Eq + neV . These scattering processes 
have the amplitudes s^o, shown in the figure. Notethat all 
quasiparticles injected from the left at energies Eo + 2neV 
and from the right at Eo + (2n -f l)eV scatter on the same 
"ladder" in energy space. 

at energy En — Em + (n — m)eV is denoted (see Fig. 

This amplitude is a function of the scattering ma- 
trix of the normal conductor and the Andreev reflection 
amphtudes \vi\ . 

To access the FCS, we make the first important ob- 
servation that quasiparticle scattering in a voltage biased 
superconducting junction is formally identical to scatter- 
ing in a normal voltage biased junction, with an applied 
harmonic ac-field. The FCS in such a system was in- 
vestigated in detail by Ivanov and Levitov[4]. Following 
Ref. [4] we note that for measurement times much longer 
than the quasiparticle scattering time, the inelastic single 
mode scattering problem can be mapped onto an elastic 
scattering problem with many modes. The scattering be- 
tween all energies and quasiparticle types of a "ladder" 
(see Fig. ^ is correlated, while different ladders con- 
tribute incoherently. We denote the ladder by its energy 
Eq in the interval —A — 2eV < Eq < —A counting its 
"leg" in the left lead (e.g. Eq in Fig^). Thus we may con- 
centrate on the scattering matrix S, with elements s"^, 
of a single ladder, and then integrate over the ladder en- 
ergy Eq. In general S has infinite dimensions, but in 
our case the vanishing probability of Andreev reflection 
far outside the gap naturally cuts the number of rele- 
vant modes to the order of A/eV. Quasiparticle current 
(but not charge current) is conserved in the scattering 
processes at the normal-superconductor interfaces. As a 
consequence S is unitary 

We are then, in line with 0, , able to directly write 
down the characteristic function in terms of all different 
many-particle scattering probabilities 

xZiA) = ^eKS"«o^"^-E,„...^"^")p^l^, (1) 

i,o 



where A is the set of counting fields A„q, one for each 
mode na. The outer sum runs over all possible sets 
of incoming modes i = {miai, TO2Q!2, . . .} and outgoing 
modes o = {ni/3i, n2/32i ■ • ■}■ 

The many-particle scattering probabilities are given by 

Aio=kri'n/(^™) n[i-/(^™')], (2) 

where f{E) = (1 -I- exp[E/kT]) is the Fermi distribution 
function and \sf\ is the determinant of the matrix formed 
by taking the columns i and rows o of S. 

Eqs. ^ and ||2J) gives us the FCS for quasiparticle 
transfer in a voltage biased superconducting junction. 
The object of main interest is however the FCS of the 
charqe transfer. We then make the second important ob- 
servation that for measuring times much longer than the 
inverse Josephson frequency, the net charge transfer is 
directly related to the quasiparticle transfer. A quasi- 
particle of type a incident at energy Em, which is scat- 
tered into an outgoing quasiparticle of type (3 at energy 
En, transports exactly m — n electrical charges across 
the junction (for m < n, the transported charge is thus 
negative for the bias in Fig. QJ. This is independent of 
the type of quasiparticles a and /3, and also of the charge 
of the quasiparticles. This can be shown by studying 
any quasiparticle scattering path (see Fig. QJ, keeping 
in mind that the process of Andreev reflection transfers 
exactly two electrons across the normal-superconductor 
interface, while a normal transmission transfers exactly 
one electron. 

This can also be seen from energy conservation: An 
electron traversing the normal part of the junction from 
left to right will absorb the energy quantum eV from 
the electric field, while an electron moving in the op- 
posite direction will emit the quantum eV. The effective 
number of quanta eV absorbed in a quasiparticle scatter- 
ing process scattered thus equals the number of electrons 
transferred from left to right. We emphasize that this 
approach correctly counts all the electrons transferred, 
including the electrons entering the superconductor as 
Cooper pairs at energies within the gap. 

We can thus count the transferred electrons by count- 
ing the transferred quasiparticles weighted by the number 
of electrons transferred in each quasiparticle scattering 
event. Thus by chosing the counting fields A„q = nX 
in Eq. we can directly write down the characteristic 
function for charqe transfer 

i,o 

Following Ref jj], the characteristic function can be 
written on a determinant form as 

XEoW=det[l^f + fSzlSx], (4) 
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where the elements of the matrix S\ are given by 
{Sx)'^i =_s^'^e*^("-™)/2 and the diagonal matrix / has 
elements = f{En)- This results provides a general 
solution to the temperature and voltage dependence of 
the full counting statistics of voltage biased supercon- 
ducting junctions, for measurement times much longer 
than inverse Josephson frequency, and the quasiparticle 
scattering time, as defined by the energy dependence of 
the scattering matrix ^ In js^'^l. Eqs. © and Q 
shows that the charge is transferred in correlated quanta 
of multiple electron charges. 

At low temperatures A;T <C A all quasiparticle states 
below the gap are filled, while all states above the gap 
are empty. This simplifies Eq. ^ significantly since it 
fixes the set of incoming modes i to all modes below the 
gap, but still all different sets of outgoing modes should 
be considered. 

The MAR ladder forms a single mode for transport in 
energy space The scattering amplitudes s^J^ can 

be decomposed into amplitudes for entering the normal 
region, t™^, propagation along the MAR ladder, t 
and leaving the normal region, tj^^, as 



T/i 



„l3a _ ^m+A ,nl3 



j-'m-t-j-l f" 
ma mn n- 



''mn — ''mk^'kr. 



{m < k < n), (5) 



and analogously for £"„ < Em- Using this decomposition 
we find that choosing two or more outgoing modes above 
the gap in the set o gives Pi|o = 0, since the matrix 
sf then contains two or more parallell rows. In other 
words, the many-particle scattering process of two or 
more quasiparticles passing the gap through the single 
mode in energy space is prohibited by the Pauli exclu- 
sion principle. For the outgoing sets with all outgoing 
modes below the gap \sf\ is evaluated through making 
sf unitary by inserting a row and a column with the 
amplitudes t^ttl^g, where Eg is the first energy above 
the gap in the ladder. Finally, the sets of one outgoing 
mode above the gap, using similar manipulations, give 
the single particle scattering probabilities. The charac- 
teristic function can then be expressed in single-particle 
scattering probabilities 



XEoW = 



m<0,ri>0 



I — m)A ^ ^ I 

a, 13 



-,0a |2 



(6) 



This simple expression is the main result of this paper. 

The cumulants of the charge transfer distribution func- 
tion are obtained by taking derivatives of hiXEoWj 
summing over all ladders (integrating over the energy 

Eo), 



{{nirD) 



-A-2eV 



-A 



dEo i-^^xrlIiXE„iX)\,^,, 

(7) 

where t is the measurement time. 

The integrands for the various moments are conve- 
niently expressed in the physically relevant n-electron 
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FIG. 2: The probability P„ = l/2eV J_^_2,y Pn{Eo)dEo of 
n-electron scattering events, as a function of voltage for n € 
{1,2,3,4} in a fully transparent junction (left) and a tunnel 
junction (right). 



scattering probabilities |l7l |. 

PniEo) = J2 QiEm+n~A) ^ \s^^^„^jEo)\' {u > 1), 

m<0 a,/3 

(8) 

with Pa(Ea), the probability of no quasiparticle passing 
the gap, giving the correct normalization J2^=o Pn^Eo) — 
1. In Fig. [3 the probability of different n-electron scat- 
tering processes, as a function of voltage, are shown. 
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FIG. 3: The zero-frequency current noise (upper panel) and 
the third cumulant (lower panel) as a function of voltage for 
transparencies D = 0.1, 0.5, 0.9, 0.99 and 1. All curves are for 
a single mode junction at temperatures kT <C A. 

The integrands for the first three cumulants in Eq. {T)) 
are the spectral current density I{E) — J2^=o ''^-Pn{E) , 
the noise spectral density Si{E) — J2^=o ''^^ ^^i^) ~ 
I{EY, and the spectral density of the third cumulant 
C^{E) = Er=o"'^"(^) - m{-iSi{E) + I{Ef). The 
expression for the current is identical to the one presented 
in Ref. jl7l |. and the expression for the zero- frequency 
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current noise reproduces the known results 0| (see the 
upper panel in Fig. |3J • 

The noise and the third cumulant, calculated numeri- 
cally, are plotted in Fig. Ofor different transparencies of 
the normal contact. We see that the subgap structure is 
generally more pronounced in the third cumulant, com- 
pared to the current and noise. Furthermore C3 changes 
sign both as a function of voltage and transparency. 

In the tunnel limit D 1 where all scattering proba- 
bilities in Eq. (jHJ are small, the corresponding probabil- 
ity distribution is Poissonian, i.e. describes an uncorre- 
lated transfer of quantas of multiple charge. Within the 
n-th subgap region, i.e. n — I < 2A/eV < n, n-electron 
transfer dominate, giving S = 2enl and C3 = {2en)^I, 
see the upper panel in Fig. 0] As seen in the right panel 
of Fig. 13 the picture is more complicated on the borders 
between the regions. Due to resonances up to three adja- 
cent n-particle processes have similar strength ^'^id 
the FCS is no longer Poissonian. 
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FIG. 4: (Upper panel) The Fano factor S/I (solid) and C3// 
(dash-dotted) in the tunnel limit D = 0.01. (Lower panel) 
The current (solid), noise (dash-dotted) and third cumulant 
(dashed) for a short, diffusive superconducting junction. 

So far we considered only single mode junctions. In 
the limit of a short junction, when the scattering ma- 
trix of the normal conductor is independent on energy 
on the scale of A, it is possible, just as for the current 
and noise to write the generating function in terms 
of the transmission eigenvalues Z3„ of the normal con- 
ductor only. For all mesoscopic conductors where the 
transmission eigenvalue distribution is known, the FCS 
can thus be obtained via averaging the single mode result 
in Eq. I^. As an example, we present in the lower panel 
of Fig. ^the first three moments for a short, diffusive 
junction. We note that the third moment is positive for 
all voltages. Just as the noise, the third moment has a 
divergency for eV 0. However, in an experiment, we 
expect the fluctuations at small voltages to be strongly 
suppressed by inelastic scattering , effectively remov- 



ing the divergency |l5| . 

In conclusion, we have derived an expression for the 
full counting statistics of a voltage biased superconduct- 
ing junction, describing a correlated transfer of quantas 
of multiple electron charge. It is found that the counting 
statistics can be expressed in terms of probabilities for 
quasiparticles to scatter between the two superconduc- 
tors. 
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